Dispersionless transport in a washboard potential 
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We study and characterize a new dynamical regime of underdamped particles in a tilted washboard 
potential. We find that for small friction in a finite range of forces the particles move essentially 
nondispersively, that is, coherently, over long intervals of time. The associated distribution of the 
particle positions moves at an essentially constant velocity and is far from Gaussian-like. This new 
regime is complementary to, and entirely different from, well-known nonlinear response and large 
dispersion regimes observed for other values of the external force. 
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Particle transport and diffusion in periodic potentials 
at finite temperatures has been addressed in so many 
contexts and over so many decades that one might think 
this to be a fully solved problem [l[. However, as modern 
experimental and numerical methods continually evolve, 
ever broader parameter regimes and time regimes become 
accessible to inquiry, and behaviors continue to be re- 
vealed that have notpreviously been explored or even 
noted [11,0,!, BBSS Ell- Intermediate time regimes 
are especially challenging. On the one hand, numerical 
methods need to be efficient to reach beyond relatively 
short time behavior. On the other, analytic methods usu- 
ally deal with asymptotia. Yet, experiments often involve 
intermediate time regimes. In this Letter we report an 
unexplored dynamical regime, namely, particle transport 
that is essentially nondispersive or coherent over long 
time intervals. 

Consider a particle moving in a periodic potential V(x) 
of amplitude Vq and period A, with coefficient of friction 
r\ at temperature T, and subject to a constant external 
force /. The variables x and t respectively denote the 
position of the particle and the time. The equation of 
motion reads 



-V'(r) - 7 r + F + ((t), 



(1) 



where r = x/\, r = (Vb/m) 1 / 2 i/A, the dot and prime 
denote derivatives with respect to r and r respectively, 
V(r) = V(x)/Vq, and the noise obeys the fluctuation- 
dissipation relation (C( T )C( r ')) = 2jTS(t — r'). Equa- 
tion fTJ) models the translational Brownian motion of a 
particle in a tilted periodic potential, and also the rota- 
tional Brownian motion of a damped pendulum driven by 
a constant torque. The pendulum provides the mathe- 
matical background underlying a number of applications, 
including mobility in superionic conductors, dynamics of 
charge-density waves, ring laser gyroscopes, and phase- 
locking phenomena in radio engineering. Perhaps most 



directly relevant for experimental testing, it also models 
a resistively and capacitively shunted single Josephson 
junction. This latter correspondence has been invoked 
in some of the most recent work on transport in tilted 
periodic potentials 0, Q. An excellent table indicat- 
ing the precise translation between the parameters of a 
number of physical systems including the quintessential 
Josephson junction testbed and our model equation can 
be found in 0- 

There are three independent parameters in the model, 
the scaled force F = Xf /Vq, the scaled temperature T = 
ksT/Vo, and the scaled dissipation 7 = -qX/ {mVo) 1 / 2 . In 
our subsequent numerical simulations we choose V(r) = 
-(1/2) cos(27rr), 7 = 0.04, and T = 0.2. 

The external force F (which we choose positive) tilts 
the potential, V(r) — Fr. For small forces, the tilted po- 
tential has wells and barriers ("washboard potential"). 
Beyond a critical force F c — tt the maxima and min- 
ima disappear, and the potential decreases monotonically 
with increasing r. For small F the motion of a particle 
is essentially confined to the potential wells from which a 
thermal fluctuation occasionally causes it to emerge. The 
velocity and dispersion of the particles in this "locked 
state" can be understood from well-known Kramers ar- 
guments. On the other hand, when F > F c , the peri- 
odic portion of the potential becomes unimportant; the 
particles follow the force with an average velocity that 
approaches v — F/j with increasing force, with a dis- 
persion given by the Einstein law of diffusion. This is 
the "running state." When the force is not too low but 
still well below F c , there is a regime where particles read- 
ily alternate between the locked and running states. As 
a result, the mean velocity of the particles rises rapidly 
with increasing force, the mobility exhibits a peak asso- 
ciated with this rapid rise, and the dispersion also ex- 
hibits a peak that reflects the bimodal distribution of 
velocities P, Q . These well-known behaviors for under- 
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FIG. 1: Mean velocity v (left axis, open circles) and mobility 
fi (right axis, black circles) vs the external force. Inset: dif- 
fusion coefficient D vs external force. The "empty" portion 
of the diffusion coefficient curve indicates that we did not 
achieve a stationary state up to the longest simulation times 
(r ~ O(10 7 )). This is the regime of persistent dispersionless 
transport (see text). 



FIG. 2: Dispersion vs time for different forces. Note that the 
curves for small (F — 0.04) and large (F = 4.0) forces quickly 
cross over from the short-time superdiffusive regime typical of 
low friction to the asymptotic diffusive regime. Intermediate 
forces exhibit a long dispersionless regime, eventually followed 
by normal dispersion. 



damped particles are illustrated in Fig. [TJ where we show 
the asymptotic mean velocity v = lim r ^ oc (r(r))/r of 
the particles as a function of F, and the associated mo- 
bility \i = dv/dF. In the inset we display the diffusion 
coefficient D = lim T ^ 00 (Ar 2 (r))/2r, where (Ar 2 (r)) — 
([r(r) — (r(r))] 2 ) is the relative mean square displacement 
or dispersion. The mobility and diffusion coefficient peak 
at a 7-dependent force more than an order of magnitude 
below F c for 7 = 0.04. Beyond the peaking force, the 
mobility reaches a constant close to I/7. The mean ve- 
locity rises quickly as particles emerge from the well, and 
settles to its steady state value in times of O(10 2 — 10 3 ) 
(explicitly estimated later) for our parameters. 

The as yet unexplored behavior occurs for forces above 
those of the mobility and diffusion peaks but below the 
critical force. In this regime, which extends from well 
below up to nearly F c , we have exhibited no values for 
D in the inset of Fig. [TJ This is because for times that 
are as long as we are able to reasonably run our simu- 
lations (t ~ 10 7 ), (Ar 2 ) has not yet reached a behavior 
proportional to time. For this wide range of intermediate 
forces (Ar 2 ) remains essentially constant over time inter- 
vals that can extend over several decades, i.e., the trans- 
port is seemingly dispersionless. This coherent behav- 
ior occurs only in the underdamped system and requires 
there to be a washboard potential. The nondispersive 
behavior is remarkable because it occurs in the presence 
of thermal fluctuations. At sufficiently long times (Ar 2 ) 
begins to grow linearly in time, presumably also for those 
cases in which we were not able to reach this regime. 

Figure[2]shows the dispersion for different values of the 
force. For very small forces (including zero force), and 
also for very large forces, (Ar 2 ) settles into the expected 
linear growth with time on the same time scale as the 



mean velocity. The corresponding diffusion coefficients 
as reflected in the linear growth slopes are those shown 
in the inset of Fig. [TJ In Fig. [2] we show two forces (F = 
2.2 and 2.5) which exhibit the full range of behaviors 
of the dispersion from early rise to dispersionless regime 
followed by normal Einstein behavior. We also show an 
example of a force, F = 0.5, for which the dispersionless 
regime extends over almost the entire time range of the 
figure after the initial rise. 

To explore the nondispersive behavior in more detail, 
in Fig [3] we display simulation results for the distribu- 
tion of the moving particles. Initially all the particles 
are placed at the minimum of one potential well, with a 
Maxwell-Boltzmann distribution of velocities at tempera- 
ture T. At early times (upper left panel) the distribution 
has a large peak around the origin because most of the 
particles are still in the initial well. As time proceeds, 
particles leave the well, giving rise to the second peak. 
The growing dispersion in this regime is due to the bi- 
modal distribution. The maximum of the second peak 
moves practically from the very beginning with velocity 
v « F/j. In time, the first peak disappears and the sec- 
ond peak settles into the behavior shown in the lower 
panels, which display the remarkable intermediate-time 
nondispersive behavior. The distribution of positions is 
a one-sided exponential that moves coherently without 
distortion or dispersion over a long time interval. This 
is particularly remarkable since during all this time (as 
well as before and after) the particles are persistently sub- 
ject to thermal fluctuations. At sufficiently long times, 
dispersion resumes and the distribution not only broad- 
ens but approaches the characteristic equilibrium Gaus- 
sian shape (upper right panel) . We have ascertained this 
same behavior for other forces. We proceed to identify 
the physical mechanism and to estimate the forces and 
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FIG. 3: Unnormalized distribution of particle positions at 
different times for F = 0.12. Top left: early times regime. 
Values of r: 500 (dot-dashed); 1000 (solid); 2000 (dotted); 
4000 (dashed). Bottom left: dispersionless time regime for 
r = 7000 (solid), 30,000 (dotted). Top right: long-time dis- 
persive regime for r = 2 x 10 5 (solid), 10 6 (dotted). Bottom 
right: Log plot of bottom left distribution to display the ex- 
ponential tail. The straight line is the exponential form @ 
with l = (F/j)t = 3600. 



crossover times for the occurrence of the dispersionless 
behavior. 

Initially all the particles are placed in one potential 
well, from which they emerge according to the well- 
known exponential distribution of first exit times r e , 
P(r e ) = Tq 1 cxp(— t c /to) [11]. The mean exit time 
To depends on the system parameters, the tempera- 
ture, and the external force. Since the particles emerge 
with velocities very narrowly distributed about v 
F/'y, this temporal distribution can be transformed to 
a spatial distribution of particles at r = 0, p(l e ,0) = 
Iq 1 exp(Z e /Zo), — oo < l e < 0, via the change of variables 
l e = (F/"f)T e . Here Iq = l To = {F/^)tq. This distribution 
looks precisely like the ones in the lower panels of Fig. [3] 
At time r the distribution will have moved to the right 
pretty much undistorted, 



P( l e,r) = ^exp [ e 



l T = (F/ 7 )r. (2) 



The moments of this distribution are consistent with the 
numerical results. In particular, when l T 3> Iq the first 
moment — > l T . The dispersion settles to the constant 
value Iq. These results reproduce the intermediate time 
behaviors highlighted in the simulation results. 

The parameter tq defines the crossover time from short 
to intermediate time behavior. Its value can be extracted 
from the simulation results. We find To = 1200 (F = 
0.12), 490 (F = 0.5), 17 (F = 2.2), 3.8 (F = 2.5). A 
theoretical estimate follows from the well-known result 
for the escape time of an underdamped particle from a 



potential well [lj, M, LLJ , r oc (T/jU)expU/T. The 
proportionality factor varies depending on the precise 
definition of To, but is generally of 0(1). For small F, 
U = 1 — F/2 + 0(F 2 ); for larger F, U is determined from 
a transcendental equation. This estimate leads to results 
consistent with the simulation values. 

What marks the existence and termination of the dis- 
persionless transport regime? There are two contribu- 
tions to the width of p(l e ,r). One is the width Iq due 
to dispersion of emergence times from the potential well, 
and the other comes from the thermal fluctuations. The 
total square width of the distribution is obtained from a 
convolution of these two contributions and is thus sim- 
ply their sum, (A?- 2 (t)) = Iq + 2Dr. We stress that D 
is the actual diffusion coefficient, which is in general dif- 
ferent from T/j (see inset in Fig. [T]). The effects of the 
thermal fluctuations are therefore submerged until times 
t ~ Iq/2D. This estimate leads to the dispersionless end- 
ing times t = 5.4 x 10 6 (F = 0.12), 8.9 x 10 4 (F = 2.2), 
and 5.6 x 10 3 (F = 2.5), consistent with the results in 
Fig. [2] The seemingly anomalous dispersionless behavior 
is thus a natural consequence of the motion of an under- 
damped particle in a washboard potential. 

A necessary condition for nondispersive motion is that 
the motion of a single particle be nonchaotic in the ab- 
sence of noise. An upper bound on the forces that lead 
to nondispersive motion is F — F c . A lower bound 
can be estimated from Fig. Q] as the force above which 
v s» F/j, which for 7 = 0.04 is around F ~ 0.1. 
To estimate this bound analytically and, more gener- 
ally, to analyze the behavior of single particle trajecto- 
ries through the dispersionless and into the dispersive 
regimes, we define the auxiliary variable z(t) that de- 
scribes the deviations of the particle trajectory from the 
approximate mean trajectory, z = r — [l r — l e ]. The 
distance l e quickly becomes negligible. The equation of 
evolution for z(t) is z = — 7rsin27r(z + l T ) — jz + C( r )i 
which describes a particle moving in a cosinusoidal po- 
tential that is itself rapidly and periodically oscillating 
in time. In the absence of noise, suppose that this par- 
ticle is near the bottom of a well, slowly relaxing toward 
the minimum. At a time "f/2F later it finds itself near 
the top of the potential, from which the small damping 
slowly pulls it toward the nearest minimum (either the 
same it was in or an adjacent one, depending on its pre- 
cise location). Before it has gone far the particle again 
finds itself near the minimum of the well as the potential 
flips, the very same minimum in which it started. Again 
the damping pulls it toward the bottom, and the cycle 
restarts. The resulting z(t) is thus expected to remain 
very small. An approximate solution for small displace- 
ments in the absence of the noise is easily obtained by 
approximating the nonlinear problem as a dichotomous 
succession of linearized problems, z = —(+)2ir 2 z — jz for 
t 2 „ < t < t 2 „+i (t 2 „+i < t < t 2 „ +2 ), where r n = nF/2j 
and n — 0, 1, 2, ■ • • . The solution of the two equations is 
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FIG. 4: Bottom panel: Two individual ^-trajectories for 
F = 0.12. The two trajectories emerge from the initial po- 
tential well at slightly different times of order 10 3 and move 
essentially nondispersively for a long time, until dispersive 
motion eventually sets in. Upper panel: associated disper- 
sion curve. 



straightforward, particularly at times ti 



z(*2n) = ae 



cos \ 9 + n arctan 



VT^A 2 
A 



, (3) 



where A = cosh u cos u, u = nj/y/2F 7 z(0) = a cos 9, and 
we have retained only leading contributions in 7/ F. The 
switching thus leads to a decay of solutions down to zero, 
which is the stabilizing behavior that leads to dispersion- 
less motion. The solution ([3]) holds when A 2 < 1, which 
places the constraint F > I.27 on the parameters. For 
the value 7 = 0.04 used in our simulations, F > 0.05, in 
gratifying agreement with the lowest forces found to lead 
to dispersionless transport. 

A clear manifestation of this description can be seen in 
Fig.QJ We show two individual z-trajectories that emerge 
from the potential well at different times. The associated 
particle trajectories are therefore in different portions of 
the exponential distribution of Fig. [3] The remarkable 
point is the flat z-trajectory seen in each case in the time 
interval between its emergence from the well to around 
t > 10 5 . This is the regime of dispersionless motion, 
as seen in the associated dispersion curve in the upper 
panel. 



These arguments do not include the effects of noise, 
which buffets z and eventually causes it to increase be- 
yond the validity of our linearized theory. The system 
then no longer finds itself near the top and bottom of 
the flip-flopping potential but rather is sufficiently out of 
phase with the extrema so that the alternation does not 
repeatedly drive it back to the same minimum. Disper- 
sion then sets in. Figure 3] confirms that the dispersion- 
less motion comes to an end, and noisy motion associated 
with linear growth of (Ar 2 ) with time resumes. 

In conclusion, we have identified and characterized a 
transport regime in the motion of thermally agitated clas- 
sical particles that involves essentially dispersionless mo- 
tion over several decades of time in appropriate parame- 
ter regimes. This remarkably coherent behavior requires 
the presence of thermal fluctuations, and is restricted 
to underdamped systems in a washboard potential. We 
have provided a theoretical framework to estimate the 
times of onset and duration of dispersionless transport 
and the parameter regimes in which it occurs. Even- 
tually, dispersive motion sets in again and the system 
resumes the familiar behavior, but it may take an inor- 
dinately long time to do so, so long that it lies beyond 
our computational means. We are confident that the re- 
sults obtained here are amenable to experimental veri- 
fication. Apart from the possibility of a mechanical or 
electronic realization of a tilted washboard system, the 
regime might be observed in Josephson junctions, where 
the voltage (the observable) is proportional to the time 
derivative of the phase (corresponding to the particle co- 
ordinate in our model). It is important to note that we 
have also observed nondispersive behavior in numerical 
simulations in two dimensions with a separable potential 
of the form V(r x ,r y ) = — (l/2)[cos(27rr x ) + cos(27rr y )]. 
This suggests that our results may be testable on the 
motion of molecules and clusters on surfaces [see refer- 
ences in [HJ]. These observations give us reason to ex- 
pect that this coherent dispersionless behavior should be 
experimentally accessible. 
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